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ABSTRACT 

The discovery of hundreds of young, bright stars within a parsec from the massive black hole Sagittarius A* at 
the center of the Galaxy presents a challenge to star formation theories. The requisite Roche densities for the 
gravitational collapse of gas clouds are most naturally achieved in accretion disks. The water maser sources 
in Keplerian rotation in the nuclei of NGC4258, NGC1068, and the Circinus galaxy indicate the presence of 
warm, extended, molecular accretion disks around black holes similar in mass to Sgr A* . Here we argue that the 
current conditions in the maser nuclei, and those near the Galactic center, represent two consecutive, recurrent 
phases in the life cycle of the nucleus of a typical gas-rich spiral bulge. The warm molecular disks that give rise 
to the observed maser emission fragment into stellar-size objects. The stellar masses, their orbital geometry, 
and the total number of stars thus formed are consistent with the values identified at the Galactic center. The 
stars tend to form in compact groups resembhng the IRS 13 complex that dominates the stellar light in the 
neighborhood of Sgr A*. 

Subject headings: accretion, accretion disks — Galaxy: center — galaxies: nuclei — masers — stars: forma- 
tion 



1. INTRODUCTION 

The dynamical sphere of influence of the massive black 
hole (MBH) at the Galactic center, measuring 1 pc in radius, is 
host to a lar ge number of luminous, massive, and h ence young 
stars (e.g., Krabbe et al. 1995; Genzel et al.'2000). The clos- 
est star to the MBH with detailed proper motion data, SO-2, 
has an apocenter distanc e of only ^ 0.01 p c and an eccen- 
tricity of ~ 0.88 (e.g., S chodel et alJl2002l) . It is a B0-O8 
dwarf with a ma ss of lO-lSM© and an age of < lOMyr 
iGhez et alJl2003at) . A number of other stars at radial dis- 
tances< 0.03 pc exhibit properties similar to those of SO-2 
(e.g., ^hez et al. 2003b). At a projected distance > 0.1 pc 
from the MBH, a cluster containing seven blue supergiants, 
known as th e IRS 13 complex, dominates the ionizing lumi- 
nosity (e.g., iNaiarro et al.lll997l) . At somewhat larger dis- 
tances of 0.1-0.5pc, as many as ~ 40 Wolf-Rayet stars 
have been observed (Genzel et al. 2003). Since the zero- 
age main sequence (ZAMS) mass of the Wolf-Rayet stars is 
30- IOOM0, this implies that the Galactic center had been 
host to a recent starburst with a total initial stellar mass of 
greater than IO^Mq. 

Various dynamical mechanisms by which young stars could 
have migrated from larger r adii into a fract ion of the central 
parse c have been considered (lGerhardl200U lGould & Ouillen 
l20()3t iKim & MorrisI l2()03l) . None are successful at ex- 
plaining the observed concentration of young stars, un- 
less an additional, hitherto undetected dynamical compo- 
nent, such a s an intermediate-mass black hole, is invoked 
dHansen & M ilosavhevic 2003). This conclusion, therefore, 
suggests that the stars have probably formed in-situ. Molecu- 
lar cloud densities of 10'"* cm~^ at the radius corresponding to 
the apocenter distance of SO-2, and those of 10''- 10'" cm""* 
at the orbital radii of the Wolf-Rayet stars, are required for 
gravitationa l collapse to o ccur in the presence of the Sagittar- 
ius A* tide ('M orrislll993h . The densities of molecular cloud 
fragments found in the central 2 pc of the Galaxy, however, 
do not exceed ^ 10'' cm"^ ('e.g.. Uackson et al.ll993l) . 



The densities necessary for gravitational collapse in the 
tidal field of the MBH could easily be sustained in accretion 
disks. The young stars may be connected to a past accre- 
tion disk at the Galactic center since their orbits appear to 
lie wi thin one or more distinct planes dLevin & Beloborodovl 
I2003t [Genzel et al. 2003). Multiple steUar disks have been 
observed in the nuclei of other galaxies at the resolution lim- 
ited radii of ~ 20 pc (e.g., Pizzella et al. 2002). 

The proposal that stars form in the outer parts of ac- 
cretion disks of active galactic nuclei (AGNs) is not new 
(Kolvkhalov & Sunvaev 1980). Th e disks are susc eptible to 
self-gravitating instability when the iToomrd (IT96? ) parame- 
ter Q = CsVL/ttGYj is less than unity, where Cj is the sound 
speed, Vl = (GMb^^/py^-^ is the angular velocity (assuming a 
Keplerian potential), Mbh is the mass of the MBH, and S is the 
surface density of gas. For a range of disk parameters, the disk 
fragments into stellar-size clumps. Shlosman & Be gelmanl 
( 1989) discussed star formation in cold, molecular disks of 
AGN. Grains in the surface layer of a disk absorb the inci- 
dent optical-to-soft X-ray flux from the AGN and reradiate it 
in the infrared. The disk bulk temperature remains low; e.g., 
if the 4 X 10''Mq black hole at the Galactic center were ac- 
creting at the Eddington rate during a past accretion episode, 
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it would have heated the disk to temperatures ~ 36r3 K, 
where r = 0.3r3 pc is the distance from the MBH. Such a cold 
disk is self-gravitating at surface densities > 5 g cm"^. The 
apparent inevitability of fragmentation implies a crisis for the 
extended disk models (.GoodmarL2003>) . 

While the construction of extended disk models may be 
challenging, nature offers evidence that these disks indeed ex- 
ist. Maser emission in the rotational level transition 6 le — > 523 
of water has been detected at 22 GHz in the nuclei of several 
Seyfert II and LINER galaxies. Keplerian rotation patterns 
and disk-like geometries of the maser sources at distances 
> 0.1 pc from the MBH have been detected in (at least) three 
cases: NGC4258 (Mivoshi et al. 1995; Greenhill et al. 1995), 
NGC1068 dGalhmore et alJll996t) . and the Circinus galaxy 
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Table 1 . 



Galaxy 


Mbh (1O'M0) 


'■disk (pc) 


o-stciiar (kms ') 


NGC4258 




0.16-0.28" 


167f 


NGC1068 


1.0* 


> 0.6" 


177^ 


Circinus 


0.17'' 


0.11-0.4'' 


167"^ 


Milky Way 


0.4"= 




75-125 



IHermstein et al.'^T999l) 
•^Greenhill et al. 1 1996) 
' Gallimore et al. (2001) 
'i Greenhill etal. 1 2003b) 

■ iGhez et al. 1 2003b) 

IPrugniel & SimierilT^^ 



JGreenhill et alJ2003albl) . This has provided an accurate mass 
determination for the central black holes, as well as the dis- 
tances to some host galaxies (e.g., Herrnstein et al. 1999). 

In Tabled we summarize the inferred values of the black 
hole mass Mbh, the radii of the masing sources from the 
MBH rjisk, and the stellar velocity dispersions <Tsteiiar of the 
host bulges. Note that in all cases the masers are located 
within the spheres of dynamical influence of the MBH rdisk < 
GMbh/cr^teiiar' which reflects the selection of these sources for 
nearly-Keplerian rotation patterns. 

The purpose of this Letter is to highlight a connection be- 
tween the warm molecular disks in AGNs and the stellar pop- 
ulations at the Galactic Center In §|2] we review the thermo- 
dynamical properties of a warm molecular disk by assuming 
that the conditions in this disk are uniformly conducive to the 
production of water maser emission. In § |3] we study the 
fragmentation of the disk and the formation of stars. 

2. CONDITIONS IN THE DISK FROM THE REQUIREMENTS FOR 
MASER ACTION 

An inverted population of the maser levels, which is a pre- 
condition for maser action, requires a minimum gas tempera- 
ture of 400 K (e.g., Maloney 2002; Watson 2002). ColUsional 
pumping of the levels requires that the density of the masing 
gas be at least 10^ cm"-'. Above 10"^ cm"-', the levels are ther- 
malized, and the maser action is quenched. We assume that 
the conditions for the operation of masers are homogeneously 
represented in the disk, and we adopt the fiducial tempera- 
ture Tgas = 400 K and density Wg^s =10'' cm"^ at a distance 
'"disk ^ 0.3 pc from a central black hole of mass Mbh ^ 1O^M0. 
We defer discussing the origin of the unexpectedly high tem- 
perature of the molecular gas to §|4] 

The assumption that the accretion disk has on average the 
same parameters as the masing regions is not essential; i.e., 
the masing regions could be localized condensations in the 
disk (§ 0}. The star formation mechanism presented in § |3l 
however, is valid even at average temperatures of less than 
400 K. 

The vertical scale height of a disk in pressure equilibrium 
is H = Cs/ ft = I]//imH,ngas, where the isothermal sound speed 
of gas with a mean molecular weight /i in units of the hy- 
drogen molecular mass is related to the temperature via c,^ = 
kTgns/ piinu^ . The vertical column density of nuclei amounts to 

E/me ~ 8.5 x lO^^^i-^/^ngT^^^rl'^M'^^^ cm"^, where «gas = 
ng X lO'' cm"^ Tgas = 400^4 K and Mbh = M? x lO^M©. 

Note that very similar column densities of 10^^ cm"^ and 
4 X 10^'* cm"^ have been inferred from hard X-ray absorption 



measures in the nucle i ofNGC106 8 (iMatt et al.l2004 and the 
Circinus galaxy (Ma tt et al.lll999h . respectively. This coinci- 
dence lends support to our characterization of the conditions 
in the disk. Detailed comparison, however, is not possible 
because the inclination and the corresponding absorbing col- 
umn densities depend sensitively on the unknown degree of 
warping in the observed disks. 

3 . FORMATION OF STARS IN A WARM MOLECULAR DISK 

Since H/r ~ 0.003(r4r3/^M7)'/^ the disk is geometri- 
cally thin, and its gravitational stability can be determined 
by evaluating the Toomre parameter Q = 2.6ng^ r^^Mj. We 
find that Q drops below unity outside the critical radius of 
~ 0.42(M7/n9)'/-' pc. The critical radii coincide, to within 
the uncertainty in n^as, with the radii of the observed maser 
disks. 

Fragmentation of a marginally unstable disk produces grav- 
itationally bound objects of the Jeans mass Mjeans = Cj / G^E ^ 

3^x~^^^nq^T^^^r^'^^'^My^MQ. The Jeans length is approxi- 
mately equal to the disk thickness, which is in turn equal to the 
Hill (Roche) radius Ru = (OT*/3Mbh)''^^'" of a newly formed 
self-gravitating object of mass (hereafter a "protostar"). 
The total disk mass Trr^j^^S within the radius of the 

masers is Mdisk 1.3 x lO'^n^^^rigT^'^rl^^M'^^^MQ. This 
mass is comparable to the conservative estimates of the 
mass of the circumnuclear disk (CND) at the Galactic center 
(Mez2er, Duschl, & Zvlka 1996 and references therein). 

The maximum number of separate protostars that can 
be produced by the Jeans instability is then A^max ^ 
A^disk/A^jeans ^ 4-300 f fi^iig^ /T4Mt , whcrc / is the fraction 
of the disk mass that undergoes fragmentation. We expect that 
the latter fraction is well below unity, e.g., f ^ 0.1%- 1%. 
The density wakes associated with the Lindblad resonances 
of the protostars that formed firs t stimulate further f ragmenta- 
tion ( Armitage & Hansen 1999i lLufkin et alJ2004h . 

Gas interior to the Hill annulus | Ar| < Ru in the disk is de- 
livered by differential rotation into the protostar' s Roche lobe 
and can accrete onto the star. The final mass of the protostar 
can thus be much larger than the Jeans mass. Equating the to- 
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tal mass in the annulus InrRnT, oc mj to the mass of the pro- 
tostar and solving for the latter yield the maximum "isola- 
tion mass" (,Lissauer..l987.) Miso ~ (4/37r)Mj;{,^/M|Jj^^ to which 
the star can grow by accreting the disk material. In gaseous 
protoplanetary disks, the growth of a planet terminates when a 
gap of width > R^ opens in the disk. In A GN disks, however, 
Miso is a large fraction of Mdisk llGoodman & T an 2003). 

If multiple Jeans patches of the disk collapse simultane- 
ously, their Hill annuli grow until they overlap. The accre- 
tion is terminated either by reaching the isolation mass or by 
the global depletion of gas from the accretion onto multiple 
protostars. 

Before the protostellar mass has reached the terminal value, 
differential rotation in the disk delivers disk material into the 
Hill sphere of the protostar at the rate m^, ^ TJ^Efi, or 

..Mio-Meyr').'/^(^)"\^) 

X ("^V'V^V r J^V'' (1) 

V400Ky V0.3pcy \WMq) 
If the protostar could accept mass at this rate, 
its growth time would be given by m*/m* ~ 

2.5 X lO''^^'lWfMy\■^^f^%''^ yr. 
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The gas arriving into the Hill sphere undergoes a shock 
near the LI and L2 La grangian points and forms two streams 
around the proto star jLubow. Seibert. & Artvmowiczlll999t 
lBateetal.ll200l . The streams circularize at a fraction of 
~ 50% of Ru from the protostar The protostar is therefore 
fed from an accretion disk of its own. We have made an im- 
plicit assumption that the entire mass entering the Hill sphere 
on horseshoe orbits is inelastically captured and remains in- 
side the sphere. 

Gas in the protostellar disk must be able to dispense 
its angular momentum if it is to contract onto a cen- 
tral protostellar core. The radial extent of the protostel- 
lar disk is large, log(/?H/2/?*) > 4, where R^, is the radius 
of a star of mass on the ZAMS. In a purely gaseous 
disk, angular momentum extraction could be achieved 
via r adial transport ("a-viscosity"; Shakura & Sunjaevj 
1973h driven by the m agnetorotational instability (MRI; 
Balbus & Hawlevl ll9 98^. via the magnetic b reakin g (MB) 
mechanism (e.g.. iMouschovias & Paleologoul lT979h . or via 
the centrifugal launching of winds from the disk surface 
(Iglandford & Pavne 1982). 

MRI and MB are expected to operate in the outermost re- 
gion of the protostellar disk if the concentration of free elec- 
trons in the gas is sufficiently large to provide for an iner- 
tial coupling with a magnetic field. This concentration de- 
pends on the detailed chemistry that we do not attempt to 
analyze here. We are thus not in the position to decide 
whether the outer edge of the protostellar disk is magneti- 
cally "dead" or "active." If it is dead, the material captured 
inside the Hill sphere accumulates in a violently unstable ring 
around the protostar. If it is active, it is, in principlem possi- 
ble that a magnetically mediated angular momentum extrac- 
tion maintains gravitational stability in the protostellar disk. 
For example, Papaloizou. Nelson. & Snellgrove (2003) found 
evidence of MB while simulating the accretion onto plane- 
tary embryos embedded within a protostellar disk with ideal 
magneto-hydrodynamics (IMHD). We proceed to check for 
gravitational stability assuming, optimistically, that IMHD is 
realized in the protostellar disk. 

If MRI is the dominant angular momentum extraction chan- 
nel, and if the accretion rate is given by equation ([0, the 
disk parameters can be evaluated in the standard fashion 
(e.g.. 'Fra nk. King. & Raine .2002) . In the optically thick, 
adiabatic l imit, we ha ve a = O.OIa-2 with a-2 ~ 0.1 - 1 
(e.g., Sanoltal] l200l . The Rosseland mean opacity of 
molecular gas above 100 K and below the opa city gap at 
^ 1000 K is in the range k 2- 10 cm^ g"' ( Pollac k et al. I 
fT994- Semenovetal.l l200l. Assuming that the outer radius 
of the disk is ~ ^h, the Toomre parameter there equals 2* ^ 
0m75fi-^ '^a°_jK°-^M°'^^r^^-^mf -^'^nf '^T^-^. (All other pa- 
rameters pertain to the parent, AGN disk medium.) The outer 
edges of an a-protostellar disk are thus gravitationally unsta- 
ble in their own right. 

If MB is the dominant angular momentum extraction chan- 
nel, the angular momentum is removed in a vertical Alfven 

crossing time of the disk, t^nb ~ '/'"Vga/^. where ~ Bunif/S 
is the specific magnetic flux threading the gas in the donor 
AGN disk, pgas is the ambient gas density (here assumed to 
equal the density of gas in the donor disk), and Bunif is the net 
magnetic field threading the disk. If Bunif is a fraction ^ of 
the equipartition field (47rpgas)'''^Cj, then (f> = 2$^fl{n/ Pgas)'"'^. 
From mass and flux conservation in the flux-freezing approxi- 
mation, we get that the column density in the protostellar disk 



is related to that in the AGN disk via = E/^. Thus, as- 
suming a uniform gas temperature, the Toomre parameters of 
the two disks are related via = £^Q. Since 1 and ^ <C 1, 
we expect 2*^1- Again, the protostellar disk is gravitation- 
ally unstable. 

We have found that the protostellar disk is suscepti- 
ble to fragmentation in its own right and separates into 
multiple clumps. Th is process h as recently been iden- 
tified in simulations JBossI 12002 : iNa kamura & Li 200^ 
'Machida. Tomisaka. & Matsu motd 12004 ). While the Hill 
sphere continues to receive gas at the rate given by equation 
O, this gas does not accrete to a single protostar but to a 
group (A'gioup > 2) of protostars sharing the same Hill sphere. 
If A^gioup ^ 3, the spatial extent of the group increases because 
of the strong two-body encounters ("dynamical evaporation") 
beyond the original Hill sphere. The evaporation takes place 
on a timescale of ~ ri~'A'gioup- The final group does not share 
the same Hill sphere even if the initial group did, and it is 
tidally stretched. 

Evidence that stars form in compact groups can be found 
at the Galactic center The IRS 13 complex is located at the 
projected distance of ~ 3" « 0.12 pc from Sgr A*and has an 
apparent diameter of less than 0.04 pc. This compact stellar 
cluster contains seven blue supergiants. Recently, additional 
candidate members of this or a similar group have been dis- 
covered (Eckart et al. 2003). Assuming an inclination of IRS 
13 relative to Sgr A*of 45° and an initial Hill radius of ~ 0.01 
pc, the spatial extent of IRS 13 would correspond to a Hill 
sphere associated with a total mass of ~ 2500Mq. This is a 
fraction of the mass of the present CND and is compatible 
with the isolation mass limit for a 3 x disk that would 

form if the CND were to start accreting toward the MBH. 

The stars forming in the disk possess initially quasi-circular 
orbits, yet several of the most bound stars, including SO-2, ex- 
hibit high eccentricities e > 0.8. While the disk torques are 
not effective at generating such high eccentricities, nearby 
encounters between stars could, in principle, scatter stars 
onto radial orbits. An encounter between stars initially be- 
longing to disjoint Hill annuli induces a velocity change 
Av ~ {Gm/Riiyl^ ^ {m/M\,byl^rVL, which is much too low 
to deflect the star by of order its own velocity, ^ ril, re- 
quired for Ae ^ 1. Significantly larger velocity kicks re- 
sult from three- or four-body gravitational slingshot interac- 
tions between binary stars sharing the same Hill sphere (mas- 
sive stars in the field generically form in compact binaries; 
e.g., ^anbeveren, de Loore, & van Rensbergen 1998). In- 
deed, the circular velocity at SO-2's apoapse, ^ 1300 km s~', 
is comparable to the maximum velocity that a 15M0 star can 
be deflected by, while avoiding direct collision with a similar 
star. 

4. DISCUSSION 

So far, we have avoided explaining how the high temper- 
ature (~ 400 K) can be maintained in the maser disk. Dis- 
sipation ofJurbulence_wdfl^ 

regard ( iNeufeld & MalonevlfT99l iDesch. WaUin. & WatsorJ 
Il998h . To demonstrate this, we search for the (unreaUstic) 
value of a that would be necessary to power the masing disk 
from within. The vertical optical depth equals r = kE ^ 

lO^^/^KtigT^^'^M'^^^rl^^ . The accretion time W = irr^^/M 
is larger by ^ - 10^ than the dynamical time. Thus, the accre- 
tion r ate M = 37ro;c?E/0 may be very diffe rent at different 
radii JGammie. Naravan. & BlandforJ 19991) and need not be 



4 



related to the bolometric luminosity of the AGN. Instead, we 
assume that ngas and Tg^s are independent variables and ob- 



1 /2 —3 /2 

tam a ^ 7 50 K~^n^^T4^M-i . This requirement is orders 
of magnitude larger than the value typically found in numeri- 
cal simulations. Therefore, the standard a-disks are not suffi- 
ciently warm to form masers. 

Next, we check whether or not the protostellar radiative out- 
put could heat the disk. Part of the residual disk gas is acces- 
sible to the ionizing radiation from the new stars. Stars with 
ZAMS masses > 8Mq ionize the am bient hydrog en within ~ 
10^ yr from their initial collapse (e.g.. lYorkell986h . and create 
H II regions within their parent disk. The characteristic region 
size is given by the Stromgren radius Rs = (3JF/47rn^aB)'''-', 
where JF is the number of ionizing photons per unit time emit- 
ted by the star and as ~ (2.6-4.5) x 10"'^ cm-^ s~' is the case - 
B recombination coefficient of hydrogen (lOsterbrockll 1 989h . 
The ionizing photon number flux for a star with an effec- 
tive temperature equals w \6Tr''kT^R\vQC~^h~^ e~'^'^°l''^* , 
where hvQ = 13.6 eV. Taking the stellar luminosity w 
(15, OOOL0)ot83, where OTg = m^/SMQ and stellar radius « 



(3.2/^o)m8°•^^ we get Rsm4x 10' 
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The H II region is much smaller than the disk radius even 
if ~ 99% of the disk gas has been depleted. The resulting 
H II distribution has a "blister" geometry typically found in 
the Orion Nebula and in other star-forming regions (SFRs). 
In SFRs, maser sources are found in the shocked, dense shell 
of the blister. The isotropic luminosities of the SFR masers, 
10"^ to ILa,, are much smaller than those of the circumnuclear 



masers, 10 to 10^ Lq. This may be the result of a longer gain 
path along coherent velocity patches in the AGN disk (where 
the gas flow is organized) rather than a larger pumping lumi- 
nosity (e.g., Deguchi & Watson 1989). 

The stars could also heat the disk gas mechanically, by driv- 
ing spiral density waves that in some circumstances lead to 
spiral shocks. Such coherent, spiral shocks in the disk have 
been proposed as the sites of maser emission ( iMaoz & McKed 
1998). This explanation, however, depends on many de- 
tails of the model that we do not attempt to investigate 
here. The most promising mechanism for heating the disk 
gas may still be illumination by the X-ray flux of the AGN 
JNeufeld. Malonev. & Congel^ll994^ . The X-rays heat a sur- 
face layer in the disk while allowing the midplane to remain 
cooler. If the average midplane temperature of the disk were, 
say, ~ 5 times lower then the Jeans mass would have been, 
A^jeans ~ O.27M0. The subsequent evolution described in §|3l 
including the rapid accretion of disk gas, the subfragmenta- 
tion into multiple stars sharing the same Hill sphere, and the 
global depletion of disk gas, would occur in the same fashion. 
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